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1. INTRODUCTION 

Consider the celebrated Korteweg-de Vries (KdV) equation 

du ^ du g du d^u _ ^ (11) 
dt ^ dx ^ dx dx^ ' 

where x and t denote the spatial and temporal variables, respectively, and ry is a nonnegative 

constant [1,2] that can be chosen as or 1. The KdV equation is used to model [3-5] the 

propagation of water waves in long, narrow, shallow channels; it also arises in other areas 

such as hydromagnetic waves in a cold plasma, ion-acoustic waves, and acoustic waves in 

harmonic crystals. 

The KdV equation is one of the most well-known and most widely analyzed nonlin- 
ear partial differential equations. It has many remarkable aspects [4,5]. For example, it 
possesses traveling wave solutions known as solitons. The numerical studies on the KdV 
equation by Zabusky and Kruskal [3] led to the discovery of multi-soliton solutions, where 
the individual solitons interact nonlinearly at close distance and then move away from 
each other without changing their shapes. In their celebrated paper [6] Gardner, Greene, 
Kruskal, and Miura showed that the initial-value problem for the KdV equation can be 
solved via the "inverse scattering transform" associated with the Schrodinger equation. 
This led to the discovery that certain nonlinear partial differential equations are "com- 
pletely integrable" : They can be solved via an inverse scattering transform and they have 
some interesting common properties such as possessing soliton solutions. Lax pairs, and 
infinitely many conserved quantities. 

A pure AT-soliton solution to the KdV equation can be written explicitly as [4,5,7,8] 

= _2^!M|iE(M), (1.2) 

where r(a;; t) is the N x N matrix whose {j, I) entry is given by [9] 

— 2KjX+8Kjt+2r)Kjt 

^ji = Sji + - T ' l<J,^<iV, (1-3) 
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with 6ji denoting the Kronecker delta, the kj are N distinct positive constants corre- 
sponding to the bound states of 0), and the Cj are N positive constants known as the 
bound-state norming constants. Pure sohton solutions are trivial in the sense that the po- 
tential u{x, 0) corresponds [3-5] to a zero reflection coeflicient in the Schrodinger equation. 
There are many important ways explicit solutions to the KdV equation may help us to 
understand nonlinearity better. For example, it is of great importance [1,10] to determine 
function spaces containing the initial data u{x,0) so that u{x,t) is globally well behaved 
(i.e. does not blow up during t G [0, +oo)) or only locally well behaved (i.e. remains 
finite during t e [0, r) and blows up at some finite time r), and explicit solutions may 
help us to understand the global or local well-posedness of initial-value problems for the 
KdV equation. Explicit solutions may contribute to the development or improvement of 
numerical methods for (1.1), and they may also be useful to check the accuracy of existing 
numerical methods. 

In this paper we present a method leading to certain explicit solutions to the KdV 
equation in the first quadrant of the a;t-plane. Let us emphasize that our aim is not to solve 
the initial-boundary- value problem for (1.1) in the quarter plane. Instead, we are interested 
in producing some explicit solutions to (1.1) in terms of truly elementary functions. We 
produce certain explicit solutions to (1.1) having the form (1.2), where T{x; t) is the matrix 
appearing in (4.1). From (4.1) it is seen that T[x;t) can be constructed explicitly by 
specifying a constant P x P matrix A, a constant P-row vector C, and a constant P- 
column vector B. We also show that such solutions can equivalently be written as in (4.9). 
In fact, it is straightforward (but tedious) to verify that the right-hand side in (4.9) is 
a (formal) solution to (1.1) no matter how A, B, and C are chosen. However, as seen 
from (4.1), an arbitrary choice for A may not guarantee the convergence of the integral in 
(4.1); even if r(x;t) obtained from (4.1) exists by choosing A appropriately, that particular 
choice for A and arbitrary choices for B and C may not assure the existence of r(a;; t)~^ 
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appearing in (4.9) for x e [0,+oo) and t e [0, r) for some r > or r = +oo. One of 
our tasks in this paper is to indicate how we may choose A, B, C in order to assure the 
existence of r(a;; t) and the positivity of its determinant for aU x e [0, +oo) and t e [0, r), 
which in turn assures the existence and weU-posedness of the solution u{x^ t) given in (4.9). 

One set of possible choices for B, C corresponds to the initial values ■u(a;,0) for 
X G [0, +oo) in such a way that u{x, 0) becomes the potential belonging to a certain class 
in the one-dimensional Schrodinger equation. For example, u(x, 0) may be viewed as a 
fragment of a real-valued and integrable potential, which has a finite first moment and 
which corresponds to a rational reflection coefficient. Other choices may be possible; e.g., 
we may further require that the corresponding one-dimensional potential u{x^ 0) vanish 
identically for a; < 0. In fact, in Section 3 we outline how A^ B, C can be explicitly 
constructed from such a potential. All such choices guarantee the existence of the integral 
in (4.1) for all x e [0, -|-oo) and each fixed t. This is because such choices, as a result of 
using (3.2), assure that each eigenvalue of A has a positive real part. Hence, the matrix 
T{x; t) defined in (4.1) exists for all x e [0, -|-oo) and each fixed t, its determinant det r(a;; t) 
is continuous in t for every x G [0, -l-cxo), and also detr(x;t) ^ 1 as a; — > -\-oo for each 
fixed t. Furthermore, for such choices it is already known [11] that the resulting tt(a;,0) is 
analytic in x and det T{x; 0) > for all x G [0, -|-oo). Consequently, we have det T{x; t) > 
for all X G [0, +oo) and t G [0, r) for some positive r. It is remarkable that, for certain 
choices of A^ B, C, we can have r = -|-oo, as we see from the example in Section 5. 

As seen from the analysis in Section 3, in case the relevant reflection coefficient is 
zero, our solution u{x,t) reduces to the pure AT-soliton solution given in (1.2)-(1.3). This 
is equivalent to choosing in (4.1) and (4.9) the matrix A as the N x N diagonal matrix 
with Kj appearing in the entry, B as the column A?"- vector having the number 1 in 
each entry, and C as the row A-vector with cj appearing in the jth entry. 

This paper is organized as follows. In section 2 we mention some of the other methods 
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to solve the KdV equation and give a brief comparison. In Section 3 we provide a physical 
motivation for the derivation of our solutions and show how they may be related to some 
scattering data. In Section 4 we show how our solutions can be obtained by solving the 
Marchenko integral equation. Finally, in Section 5 we present an example to illustrate our 
method. 

2. SOME OTHER METHODS FOR THE KdV EQUATION 

As seen from Section 3, our method is based on using the inverse scattering transform, 
exploiting the degeneracy of the kernel of the Marchenko integral equation as indicated in 
(4.7), and solving the Marchenko equation (4.4) algebraically. There are also methods to 
solve the KdV equation without using the inverse scattering transform. For example, the 
technique [12,13] based on using the Backlund transformation, the technique [8,14] using 
the Darboux-Crum transformation, the Wronskian techniques and their generalizations 
[7,15-18], and the Hirota method [19]. Such methods are also used to produce certain exact 
solutions to the KdV equation. The idea behind the methods using the transformations 
of Backlund and Darboiix-Crum is to obtain new solutions to the KdV equation from 
other previously known solutions. The basic idea behind the Wronskian methods and the 
Hirota method is to represent the solution to the KdV equation in a particular form so 
that certain functions in the representation satisfy certain linear differential equations even 
though the solution itself satisfies a nonlinear differential equation. The explicit solutions 
produced by our method have the same representation (1.2) or (4.2) as in the Wronskian 
methods; however, our matrix T{x;t) (or a part of it) does not necessarily satisfy a linear 
partial differential equation as expected in the Wronskian methods. In the method based 
on the Darboux-Crum transformation, the solution to the KdV equation has the same 
representation as in (1.2) or (4.2), provided the initial solution is chosen as zero; there is 
certainly some connection between that method and our method because they both yield 
the AT-soliton solution in the easiest case; however, any possible connection in the more 
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general case is not clear at the moment and requires a detailed analysis, which we plan to 
do in the future. For the time being, we only emphasize that our exact solutions satisfy 
the half-line KdV equation with the drift term rjUx where we can choose rj = or rj > 
at will, they include some global-in-time solutions as well as some local-in-time solutions, 
and they are algebraic combinations of truly elementary functions. One advantage of our 
method is that it can be generalized to obtain certain explicit solutions to the matrix KdV 
equation as well as to the scalar and matrix nonlinear Schrodinger equations. 

Some other explicit solutions to the KdV equation known in the literature include 
algebraic solitons [20-22], rational solutions [22,23], various singular solutions [24-26] such 
as positons and negatons, solutions [22] to the periodic and other KdV equations, solutions 
[27] that are not quite as explicit but expressed in terms of certain projection operators, 
and various other solutions [28,29]. It is already known that some rational solutions can 
be obtained by letting the bound-state energies go to zero in the A'^-soliton solutions. We 
plan to do in the future a detailed comparison between our solutions (and their possible 
generalizations) and exact solutions obtained by other methods. Some generalizations of 
our solutions might be obtained by letting the dimension of the matrix A in (3.3) go to 
infinity, by choosing the entries of A and C given in (3.3) in some particular way or by 
letting some entries go to certain limits such as zero, and by analyzing the singularities 
encountered at t — t. 

3. SOME POSSIBLE CHOICES FOR A, B, C 

In this section we indicate a possible set of choices for A, B, C appearing in r(a;; t) of 
(4.1) so that the resulting function u{x,t) given in (1.2) or (4.2), or equivalently that in 
(4.9), is an explicit solution to (1.1) for all x e [0, -|-oo) and t e [0, r) with some positive r. 

Starting with the initial value u{x,0) with a; > 0, we extend it to the whole line by 
choosing u{x,0) = for a; < and we uniquely determine the corresponding scattering 
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data {R,{K,j},{cj}}. Here, R{k) is the corresponding right reflection coefficient [30-34], 
the set of constants Kj with < ki < ■ ■ ■ < kn corresponds to the bound states associated 
with the fuU-hne potential u{x,0), and the set of constants cj corresponds to the associated 
bound-state norming constants. The construction of {R, {nj}, {cj}} can be accomphshed 
through the foUowing steps: 

(a) Given u{x,0) for x e [0,+oo), uniquely determine the corresponding Jost solution 
fr{k, x) from the right by solving the initial- value problem for the half-line Schrodinger 
equation 

^ + k'f. = uix,0)fr; fr{k,0) = l, ^^^ = -ik. 

(b) Recover the corresponding right reflection coefficient R and the transmission coefficient 
T with the help of the asymptotics [30-34] of /r as a; — > -|-oo, namely by using 

Mk, ^) = ^ ^"'^ + Y§j '^'^ + ^^1)' ^ ^ 

It is known [30-34] that T is related to R via 



where C+ := U R, is the upper half complex plane, and the 0+ indicates that 
the limit from C+ should be used to evaluate T{k) for real k values. 

(c) Construct the set {kjI^Li by using (3.1). 

(d) Construct the set of positive constants {cj}f^i by using [30] 



Cj = -[Res{T,iKj)f 



1 

h / dxfr{iKj,xy 

-K'j Jo 



where the purely imaginary constant Res(T, iKj) denotes the residue of T at A; = ikj. 
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Having constructed R which is a rational function of A;, we determine all its poles in 
C"*" and the coefficients in the partial fraction expansion of R at such poles. It is known 
[30-34] that R{—k*) — R{k)* with the asterisk denoting complex conjugation, and hence 
such poles are either located on the positive imaginary axis I""" or they occur in pairs 
symmetrically located with respect to 1+. Let us use M to denote the number of poles 
in C+ without counting the multiplicities, and let us order them in such a way that the 
first n pairs are located off 1+ at A; = ±aj + if3j with aj > and < /?i < • • • < in 
case several distinct aj values correspond to the same Pj, we can further arrange aj in 
increasing order. We choose our notation so that the remaining M — 2n poles occur at 
k = iujj on 1+ with < a;2n+i < • • • < <^m- We let rrij indicate the multiplicity of the jth 
pole. 



Let Hi? denote the part of the partial fraction expansion of R containing only the 
poles in C"*". We have 



n rnj 



i=l s=l 



(k — if3j — ajY {k — iPj + Oij)^ 



i=2n+l s 



As a result of R{—k*) = R{ky, the constants e^s, 7^5, and rjs appearing in (3.2) are all 
real; in fact, we have 



{rrij — s) \ dk'^i 



[R{k) {k-aj- i(3j 



[R{k) {k-iujjY 



For j = 1, . . . , n, let us define Cj := 2 [^jm^ ^jmj 



j = l,...,n, 

j = 2n+l,...,M. 

. 7ji eji ] and 
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where I2 denotes the 2x2 unit matrix, each column vector Bj has 2mj components, each 
Aj has size 2m j x 2m and each 2x2 matrix is defined as 











A. 




—n/ ■ R- 
a-, Pj _ 




Similarly, for j 


= 2n + l,.. 


. , M, let 










-iVj 


-1 





... 


- 











iVj 


-1 


... 







"0" 












... 









Ar.= 












, Bj := 
















... iVj 


-1 




.1. 




. 








... 









where each column vector Bj has m,- components and each Aj has size mj x mj. Note that 



we can write (3.2) as 



UR{k) = -i[Ci 



M 



{k-iAi)~^ 

{k-iA2)-^ ... 








lBmj 



... {k-iAM)-^^ 

The above expression corresponds to a minimal realization [35] of IIR. Associated with the 
bound-state data {lij, Cj}jLi, we let 



A 



M+j 



Let us also define 



A := 



Ai 
A2 



Kj, C 








M+j '■— Cj, Bm+j 1, 



j = l,...,iV. 



, B:= 



L Bm+N 



, C :— [Ci ... Cm+n]- (3.3) 



... Am+nJ 
Note that ^ is a P x P block square matrix, S is a column P-vector, and C is a row 
P- vector, where P is the constant given by 

n M 

P:=N + 2^mj+ ^ mj. 

j=l j=2n+l 
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We also note that all the entries in A, B, and C are real constants. 

4. EXPLICIT SOLUTIONS 

In this section we construct our explicit solutions in terms of the three matrices A, 
B, and C. In Section 3 we have described how A, B, C may be related to some scattering 
data. Let us define 

poo 

r(x; t) ■.^Ip+ / dz e-^^SCe-^^e«*^'+'^^*, (4.1) 

J X 

where Ip is the P x P unit matrix. Our main result is that the quantity u{x, t) given as 

d 



^ detT(x;t) 
ox 

det r{x;t) 



(4.2) 



is a solution to (1.1) as long as det r(x; t) > or, equivalently, as long as the matrix r(a;; t) 
is invertible. It is known [11] that detr(a;; 0) > for a; e [0,+oo). As seen from (4.1), the 
matrix r(a;; t) can be explicitly constructed from A, B, and C, and as argued in Section 1 
we have detr{x;t) > for all x e [0, +oo) and t e [0, r) for some r > 0. There are two 
possibilities: If r = +oo then the solution u{x, t) given in (4.2) is a global-in-time solution 
to (1.1); otherwise, it is a local- in-time solution. 

The proof that (4.2) satisfies (1.1) when T{x; t) is invertible can be outlined as follows. 
The solution to (1.1) via the inverse scattering transform is obtained as in the diagram 
{i?(A;),{«,},{c,}} «(x,0) 

direct scattering 

time evolution solution to KdV ('^•^) 

{i?(A;)e«^'='*-2^^'=*, {«,•}, {c, e««^+2^«^-^}} ^"^^'^^^ u{x,t) 
The inverse scattering step in (4.3) for a; > can be accomplished by solving the time- 
evolved Marchenko equation [4,5,9] 

/"OO 

K{x,y;t) + ^{x + y;t) + dz K{x, z;t) ^{y + z;t) = 0, y > x > 0, (4.4) 

J X 
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where the Marchenko kernel ^{y; t) is given by 

-, poo N 

n(y; t):=— dk R{k) e^'^^^-'^'^i^^+'^y + y cj e8«^+2^«^*-«^2'. (4.5) 
i_oo ^ 

If t = in (4.5) then we can exphcitly evaluate il(j/; 0) in terms of A, B, C given in (3.3), 
and this can be accomplished with the help of the generalized Cauchy integral formula 
by using a contour integration along the boundary of C"*". In general, we cannot evaluate 
f2(j/; t) the same way for all t > 0, although there are cases when we might be able do this; 
for example, if all the eigenvalues of + 27] A have nonpositive real parts, then we might 
explicitly evaluate ri(j/; t) and obtain 

n{y- 1) = Ce^'^'+^'^'^'-y'^B. (4.6) 

It turns out that the evaluation of (4.5) as (4.6) yields (4.9), which is a solution to (1.1) 
as long as T{x\t) is invertible. As discussed in Section 1, this invertibility holds for all 
X G [0, +cxd) and t G [0, r) for some r > 0, where the value of r depends on the value of t] 
and the entries of the constant matrices A and C given in (3.3). We can write VL{x + y; t) 
as a dot product of a P-vector not containing x and a P-vector not containing y. This 
separability is easily seen from (4.6) by writing 

^{x + y;t) = Ce8M3+2^At-xAg-j/^^^ ^^j^ 

where Ce^*'^^'^'^'^^*~^'^ is a row P-vector and e~y^B is a column P-vector. The degeneracy 
of the kernel Q{y;t) allows us to solve (4.4) explicitly by algebraic means. In fact, its explicit 
solution is given by 

K{x, y; t) = -Ce^*^'+2''^*-^^r(a;; t)-'e-^^S, (4.8) 

where r{x;t) is the matrix in (4.1). Finally, the time-evolved potential u{x,t), which is 
also a solution to (1.1), is obtained from (4.8) via [4,5,9] 

dK{x,x;t) 
u{x,t) = -2^^, 
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leading to 



u(x,t) = 2^ [Ce«*^'+2^^*-^^r(a;;t)-^e-^^S 
ox I 



From (4.1) and (4.9) we obtain 



d 



u(x, t) = —2— tr 
ox 



T{x;t)-^^r{x;t) 



(4.9) 



(4.10) 



where we have used the fact that in evaluating the trace of a product of two matrices, we 
can change the order in the product. Using Theorem 7.3 on p. 38 of [36], we can write 
(4.10) also as (1.2) or (4.2). 

As indicated in Section 1, it is somehow surprising that any set of arbitrary choices 
for A, S, C in (4.1) and (4.9) yields a formal solution to (1.1). It can independently and 
directly be verified that u{x, t) given in (4.9) is a solution to (1.1) in a region in the xt-plane 
as long as T{x; t) exists and is invertible in that region. The verification of this can be 
achieved in a straightforward way by taking the appropriate derivatives of the right-hand 
side of (4.9) and substituting them in the left-hand side in (1.1). 

5. AN EXAMPLE 



We will now illustrate our method by an explicit example. Consider the scattering 
data with no bound states and 

-2ie{k - i/2) - V^-f 



where e and 7 are some positive constants. Using (5.1) in (3.3) we obtain 



1/2 -^3/2 
^/3/2 1/2 



B 



C = 2[7 e] 



(5.1) 



(5.2) 



Alternatively, we can start with A, B, C given in (5.2) without even knowing that they 
may be related to some scattering data. The use of (5.2) in (4.1) results in 

det T{x; t)=l-^ (e^ + 7^) e2(^-8)*-2a^ 
1 



+ e(v-8)t-x _ sin(^^i _ + (e + VSl) cos(v^r/t - y/Sx) 



(5.3) 
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Note that detT{x;t) > for all x,t>Oif (e^ + 7^) < 4/9 and < < 8. It can directly 
be verified that u{x,t) obtained as in (4.2) with detr(x;t) given in (5.3) solves (1.1) and 
hence it is a global-in-time solution. Not imposing such restrictions on e, 7, and r/, we still 
obtain solutions to (1.1), which may however be only locally well behaved or may even 
have singularities. 

For example, by choosing e = 7 = 1/2 and r] = 1, we obtain the explicit solution to 
(1.1) in the form 

(j){x,t) 



u{x, t) 



1 _ ^e-2(^+^*) + ^e-(^+^*) cos(\/3(x -t)+ 7r/12) 



2 ' 



where we have defined 

<j){x,t) := 6e-2(^+^^)-4\/2e-("+^*) sin(V3(a;-t)-7r/12)-4=e"^^^+^*^ sin(V3(a;-t)+7r/4). 

v2 

This solution is valid for all x G [0, -|-oo) and t G [0, -l-oo), and its Mathematica animation 
is available [37]. 





"0" 


, B = 


1 




1 



C=[27 2e ci]. 



(5.4) 



Adding bound states in our example results in global-in-time solutions containing 
solitons. For example, by choosing 

" 1/2 -a/3/2 
A= ^/2 1/2 

Ki 

we get another explicit solution to (1.1) valid for all x G [0, -|-oo) and t G [0,-|-oo). An 
explicit display of u{x^t) corresponding to (5.4) is available in a Mathematica file, but 
it takes many pages to display it; its animation with e = 1/2, 7 = 1/2, = 1, ki = 2, 
and ci = 3 is also available in the same Mathematica file [37] . The explicit global-in-time 
solution and its Mathematica animation are also available [37] in a Mathematica file for 
the choices 

1/2 -73/ 



A = 





L 











- 




-0- 








, B = 


1 


1 


Ki 












.1. 
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C=[27 2e ci C2], 



which contains two sohtons. 
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